SECURITY  CLASSIFICATION  OF  THIS  pace  flWion  Dora  Entormd) 


REPORT  DOCUMENTATION  PAGE 


READ  INSTRUCTIONS 
BEFORE  COMPLETING  FORM 


I.  REPORT  NUMBER 


tm 


.  pw-coaa- 

L-E  (and  Submit) 


12.  GOVT  ACCESSION  NO 


3.  RECIPIENT'S  CATALOG  NUMBER 


Technical  Memo 
Study ‘ 


19 


'  BRAGG-CELL  Receiver 


5.  TYPE  OF  REPORT  6  PERIOD  COVERED 


FINAL 


6.  PERFORMING  ORG.  REPORT  NUMBER 


D-43Q7, 


7.  AUTHORS A) 

SigPro  Systems  Inc. 
(Dr.  Lonnie  Wilson) 


S.  CONTRACT  OR  GRANT  NUMBERfA) 


NAS7 -  9  16 


».  PERFORMING  ORGANIZATION  NAME  AND  ADDRESS 

Jet  Propulsion  Laboratory.  ATTN:  171-209 
California  Institute  of  Technology 

QaJt— Ozloilb^  RASAdena.  .CA...9L1Q3 - 


10.  PROGRAM  ELEMENT.  PROJECT,  TASK 
AREA  a  WORK  UNIT  NUMBERS 


RE  182  AMEND  >187 


•  I.  CONTROLLING  OFFICE  NAME  AND  ADDRESS 

Commander ,  USAICS 
ATTN:  ATSI-CD-SF 

Hiij> nhucft  A7.  I  3 - 7 OOP 


12.  REPORT  DATE 

_ 10  Apr  87 


{4.  MONITORING  AGENCY  NAME  t  ADORESSflf  dlltarant  from  Controlling  Otllco) 

Commander ,  USAICS 
ATTN:  ATSI-CD-SF 
ft.  HuAPhura.  AZ_SS61 3--7Q.Q0- 


13.  NUMBER  OF  PAGES 
81 


IS.  SECURITY  CLASS,  (of  t hit  roport) 


UNCLASSIFIED 


IS*.  DECLASSIFICATION/ DOWNGRADING 
SCHEDULE 

NONE 


16.  DISTRIBUTION  STATEMENT  (ol  thla  Roport) 


Approved  for  Public  Dissemination 


17.  DISTRIBUTION  STATEMENT  (ol  Ola  abottdc t  an torod  In  Block  20,  If  dtdarant  Irani  Roport) 


ta.  supplementary  notes 

Frepared  by  Jet  Propulsion  Laboratory  for  the  US  Army  Intelli 
gence  Center  and  School's  Combat  Developer’s  Support  Facility. 


It.  KEY  WORDS  (Continue  on  rovroo  •<«*•  II  noc •••ary  and  identity  by  block  numbor) 

Bragg-Cell  Radar,  Frequency  Response,  IF  Signal,  Nonlinear  Acoustic 
Signal,  Computer  Simulation,  Uniform  and  Gaussian  Statistics,  AO 
Transducer 


20.  ABSTRACT  rContfaua  aa  rararaa  atda  ft  non ooaary  and  Idorutlr  by  block  numbor) 

A  first-order  mathematical  model  of  the  frequency  measuring  portior 
of  the  Bragg-Cell  radar  receiver  is  developed  in  this  report.  This 
model  is  implemented  in  a  computer  program  (listed  in  the  appendix) 
Various  combinations  of  assumed  uniform  and  Gaussian  distributions 
for  the  IF  and  acoustic  signals  are  input  to  obtain  a  graphical^ 
frequency  output  from  the  model.  Comparison  with.  a.  Ei  agg-utx 

output  is  not  attempted. 


7057-16 


U.S.  ARMY  INTELLIGENCE  CENTER  AND  SCHOOL 
Software  Analysis  and  Management  System 


BRAGG-CELL  RECEIVER  STUDY 


EAAF 


Technical  Memorandum  No.  19 


10  April  1987 


Author;-^ 

Qv  IjuaJ^ptn 


Lonnie  A.  Wilson,  Consultant 
SigPro  Systems  Incorporated 


Approval : 


es  W.  Gillis,  Subgroup  Leader 
Algorithm  Analysis  Subgroup 


dward  J.  Records,  Supervisor 
USAMS  Task 


A.'  FY  E liman.  Manager 
Ground  Data  Systems  Section 


Fred  vote ,  Manager 
Advanced  Tactical  Systems 


JET  PROPULSION  LABORATORY 
California  Institute  of  Technology 
Pasadena,  California 


JPL  D-4307 


K:v 

KW' 


kt’iVi 

m 

m 

ten; 

I 

m 


mi 

m 


m 
% 
f 


EK'W 

KWA 

i 

Vi  ,» 


tW 

•.VlSv 


$$«*• 

W|*® 

n.nli 


K# 

'His: 


m 


'.y5‘ 


m 

SviWv 

kS;w 

Bgfe 

W» 


„’»*»  w 


&s$: 

wlv'i 

A*:-: 


(SsS: 


B® 


EXECUTIVE  SUMMARY 


This  Technical  Memorandum  was  prepared 
originally  as  part  of  the  Generic  ELINT/COMINT  Sensor 
Report  (FY-86)  which  was  elmininated  under  the  FY-87 
statement  of  work  (SOW  #2)  ,  undated  (delivered  to 
JPL  19  November  1986) . 

The  purpose  of  the  Generic  ELINT/COMINT  Sensor 
Report,  of  which  this  paper  was  intended  (in  its 
final  form)  to  become  part  of,  was  to  establish  a 
basic  superhetrodyne  receiver  based  sensor  model  and 
perform  simulations  with  it  to  determine  the  shaping 
or  coloring  of  the  statistical  distributions  of  the 
radar  free-space  signal  parametrics  by  a  typical 
sensor  prior  to  reaching  the  self-correlation 
processes.  It  was  also  intended  for  incorporation 
into  the  algorithm  test  bed  so  algorithms  could  be 
tested  with  realistic  distorted  data  rather  than 
unrealistic  stastically  pure  data. 

This  work  was  originated  in  support  of 
unanswered  questions  from  previous  self-correlation 
studies.  The  modeling  and  simulation  approach  was 
used  because  "live  date"  could  not  be  obtained. 

This  paper  is  being  published  because  it  was 
completed  in  FY-86  with  FY-86  funds  and  still  serves 
a  useful  function. 

Several  significant  results  were  noted  while 
performing  simulations  on  this  model.  The  first 
order  math  model  is  developed  for  the  Bragg-cell 
receiver.  Whether  the  input  distributions  were 
uniform  or  Gaussian,  the  output  distributions  were 
found  to  be  shaped  so  as  to  append  a  staircased 
triangular  distribution. 
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ABSTRACT 


J 

Bragg-cell  receivers  are  employed  in  specialized  Electronic  Warfare 
(EW)  applications  for  the  measurement  of  frequency.  Bragg-cell  receiver 
characteristics  are  fully  characterized  for  simple  RF  emitter  conditions, 
but  less  understood  for  complex  and  wideband  RF  emitter  signals.  This 
receiver  is  early  in  its  development  cycle  when  compared  to  the  IFM  receiver. 

Functional  mathematical  models  are  derived  and  presented  in  this 
report  for  the  Bragg-cell  receiver.  Theoretical  analysis  is  presented  and 
digital  computer  signal  processing  results  are  presented  for  the  Bragg-cell 
receiver.  Probability  density  function  analysis  are  performed  for  output 
frequency. 

Probability  density  function  distributions  are  observed  to  depart  from 
assumed  distributions  for  wideband  and  complex  RF  signals.  This  analysis 
is  significant  for  high  resolution  and  fine  grain  EW  Bragg-cell  receiver 


systems. 
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INTRODUCTION 


The  development  effort  is  for  theoretical  analysis,  software  algorithms 
implementations,  signal  processing  analysis,  and  test  analysis  for  probability 
density  function  characterization  of  a  Bragg-cell  receiver  system.  SigPro 
Systems  Inc.  is  performing  this  work  under  JPL  Contract  No.  957474,  Mod. 
No.  1. 

This  final  report  describes  the  completed  study  on  the  Bragg-cell 
receiver.  A  functional  mathematical  model  is  derived  for  the  receiver. 
This  math  model  is  implemented  into  digital  computer  software  and  signal 
processing  analyses  are  performed. 

Theoretical  analysis  is  developed  for  the  Bragg-cell  receiver.  Ideal 
and  practiced  subsystem  implementations  are  presented  and  analyzed.  The 
key  ELINT  parameter  of  interest  is  frequency,  hence,  the  theoretical  analysis 
presents  the  basic  analysis  approach  for  frequency.  Probability  density 
function  analysis  is  presented  for  first-order  subsystem  modeL 

The  theoretical  model  for  the  Bragg-cell  is  developed  into  a  digital 
computer  program.  All  computer  programs  are  listed  in  Appendix  A. 

Signal  processing  analysis  is  performed  using  the  derived  first-order 
model.  The  pdf  characteristics  of  the  Bragg-cell  receiver's  output  frequency 
are  determined  for  assumed  input  frequency  pdf  (constant  and  Gaussian) 
characteristics. 


BRAGG  CELL  RECEIVER 

INTRODUCTION 

A  Bragg— :ell  receiver  is  used  to  measure  frequencies  of  input  RF  radar 
pulses,  with  wide  frequency  coverage,  large  dynamic  ranges  and  simultaneous 
or  time  coincident  signals.  The  ideal  Bragg-cell  receiver  will  measure  the 
frequency  parameter  on  a  single  pulse  basis  with  no  distortion.  Practical 
Bragg-cell  receivers  will  distort  the  frequency  characteristic  of  a  radar  emitter 
because  the  ideal  Bragg-cell  receiver  characteristics  are  difficult  to  approximate 
in  real  life. 

The  radar's  frequency  parameter  can  be  characterized  using  a  random 
variable.  The  random  variable  may  be  quantified  by  the  probability  density 
function  (pdf),  mean  value,  variance  and  other  moments.  As  a  radar  signal 
passes  through  the  Bragg-cell  receiver  the  random  signal  description  is  modified 
or  distorted  by  the  receiver's  nonlinear  transfer  function  characteristic. 
Bragg-cell  receiver  characterization  of  RF  frequency  will  differ  from  actual 
ra .  .  emitter  characterization  by  multifaceted  distortion  effects. 

This  section  presents  u  functional  description  of  the  Bragg-cell  receiver 
with  potential  distortion  sources  and  mechanisms  and  a  first-order  Bragg-cell 
receiver  model  is  developed.  This  Bragg-cell  receiver  model  Is  used  to  analyze 
the  distortion  of  the  frequency's  pdf.  The  first-order  Bragg-cell  receiver  model 
is  limited  to  key  distortion  effects. 

IDEAL  BRAGG-CELL  RECEIVER  CHARACTERIZATION 

A  functional  block  diagram  of  a  Bragg-cell  receiver  is  shown  in  Figure 
1.  The  basic  receiver  includes  an  RF-to-IF  downconverter,  Bragg-cell  with 
AO  transducer,  laser  source,  beam  expander  optics,  beam  focusing  optics, 
detector  array  and  processing  electronics. 
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RP-TO-IP 


The  RF  input  signal  is  characterized  bv 


=  XCfc)  cosfwct,  +  <j>c*5] 


E q.  1 


where 


x(t)  = 

Pulsed  or  CW  signal, 

X&> 

amplitude  modulation  function. 

UJC  - 

fixed  angular  carrier  frequency, 

phase  deviation  function. 

The  RF-to-IF  downconverter  is  described  in  Reference  1.  The  ideal 
downconverter  is  composed  of  an  ideal  local  oscillator  which  is  a  perfect  CW 
signal,  a  perfect  product  mixer,  and  an  ideal  bandpass  filter.  The  ideal  IF 
signal  output  from  the  downconverter  is  described  quantitatively  by: 


sz  x&)  cos  «+* 


Eq.  2 


Equation  2  is  identical  to  Equation  1  with  the  fixed  angular  frequency  being 
changed  from  wc  to  wI?.  The  IF  center  frequency  is  chosen  to  match  the 
center  frequency  of  the  Bragg-cell  and  AO  transducer.  Ideally  the  IF  signal 
replicates  the  input  RF  signal.  The  RF-to-IF  downconverter  is  described  and 
characterized  in  Reference  1.  RF-to-IF  downconverter  analysis  will  not  be 
presented  in  this  report. 

The  IF  signal  drives  the  AO  transducer  which  launches  an  acoustic  signal 
through  the  Bragg-cell.  Specifically,  the  IF  electronic  signal  is  converted 
into  an  acoustic  wave  (slow  wave  replication  of  the  IF  signal)  which  propagates 


'/wvww'Wi.VLN  n  y.'rMV'.’Vtftuiiywwrtv  T^iwiwwa'TOwy. 


through  the  optically  transparent  Bragg-cell.  Througn  the  eiasto-ootic  effect, 
the  acoustic  wave  produces  a  spatial  modulation  of  the  refractive  index  in 
the  Bragg-cell  (Reference  2.) 

The  laser  source  is  expanded  by  the  beam  expander  optics  and  illuminates 
the  3ragg-ceiL  As  the  laser  Lignt  passes  through  the  Bragg-cell.  the  refractive 
index  variations  produced  by  the  acoustics  wave  (which  is  a  slow  form  replica 
of  the  IF  signal)  are  impressed  onto  the  optical  signal  as  a  spatial  phase 
modulation. 

Focusing  optics  are  used  to  equivalently  Fourier  transform  the  modulated 
optical  beam.  Resultant  optical  Fourier  transform  signal  is  focused  onto  the 
detector  array  or  frequency  focal  plane  (see  Figure  1.)  The  photodiode  detector 
array  is  used  to  capture  and  convert  the  optical  signal  to  an  electronic  signal 
representation.  The  detector  array  is  a  linear  photodiode  array  which  is 
functionally  used  to  detect  and  measure  light  intensity  versus  frequency  of 
the  transformed  optical  signal.  The  detector  array  is  electronically  read  by 
the  processing  electronics,  which  also  provides  follow-on  signal  processing 
functions. 

The  Bragg-cell  modulator  is  depicted  in  Figure  2.  The  incident  laser 
light  interacts  with  the  acoustic  signal  in  the  cell  to  produce  a  diffracted  light 
information  signal  along  with  undeflected  laser  light  signal.  The  diffracted 
laser  signal  is  also  frequency  shifted  by  the  acoustic  signal  frequency.  The 
incident  laser  light  hits  the  Bragg-cell  at  an  angle  9^  with  respect  to  the 
s  axis.  The  IF  electronic  signal  is  sent  to  the  AO  transducer  to  convert  the 
IF  signal  to  a  replica  acoustic  signal.  The  acoustic  wave  of  frequency  1 1 
propagates  with  velocity  along  the  x  -  axis  in  the  Sragg^celL 

The  diffracted  laser  light  signal  results  from  the  AO  interaction  effects 
between  the  incident  laser  beam  and  the  acoustic  information  bearing  signal 
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in  the  Bragg-celL  Maximum  diffracted  laser  light  energy  occurs  in  the  first 
order  beam  when  (Reference  2): 


Eq.  3. 


<2 t\J, 


©*  =  ®* 


Eq.  4 


where 


=  Wavelength  of  laser  light 

=  Frequency  of  acoustic  wave 

=  Velocity  of  acoustic  wave  in  Bragg-cell  along  the 
x-axis. 

©  j  =  Effective  angle  of  incident  laser  light. 

*  Effective  angle  of  first  order  diffracted  laser  light 
information  signal. 

Detailed  Bragg-cell  analysis  would  require  AO  cell  refractive  index 
variations  analysis  for  precise  quantification.  This  first  order  Bragg-cell  model 
does  not  require  this  detailed  analysis.  Functionally  the  Bragg-cell  model 
is  accurate  and  conceptually  clear  for  the  external  angles  9  ;  and  9^. 

Bragg-cell  processors  are,  in  general,  configured  so  that  9.«  0.1 
radian,  hence,  Equation  3  can  be  readily  and  accurately  simplified  to: 
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Using  Equation  4  and  5,  the  following  equation  for  the  effective  angle 
of  the  first  order  diffracted  laser  light  information  signal  is: 


2VS 


Eq.  6 


Equation  6  reveals  that  9^  varies  directly  with  the  acoustic  wave  frequency. 
As  stated  earlier,  the  EF  signal  frequency  is  related  to  the  acoustic  wave 
frequency  by:  (Reference  2) 


Eq.  7. 


where  K  is  a  constant. 

Equation  8  results  by  combining  Equations  6  and  7, 


On  — 


2,  K  Vs 


Eq.  8. 


Equation  8  shows  that  9 .  varies  directly  with  the  IF  signal  carrier  frequency. 

d 

The  actual  displacement  of  the  diffracted  laser  Information  signal  from 
the  undeflected  laser  light  (see  Figure  2)  is  given  by  (Reference  2.) 


i  wtviTwisWiVii 


where  F  is  the  focal  length  of  the  focusing  optics. 

A  fixed  Bragg-cell  receiver  will  result  in  Equation  9  reducing  to: 

A=  K,  -¥-xp  mi 

The  actual  displacement,  d,  of  the  diffracted  laser  beam  relative  to 
the  undeflected  laser  light  passing  through  the  Bragg-cell  is  directly  proportional 
to  the  IF  frequency  of  the  input  IF  electronic  signal.  The  diffracted  laser 
beam  will  hit  the  linear  detector  array  at  a  specific  location  which  corresponds 
to  one  IF  frequency.  The  Bragg-cell  receiver  is  calibrated  by  varying  the  IF 
frequency  of  the  input  IF  signal  over  the  entire  operating  bandwidth  and 
capturing  and  analyzing  the  electrical  output  signal  from  the  linear  detector 
array.  Detector  array  element  signal  versus  IF  frequency  characteristic  is 
used  to  quantify  the  Bragg-cell  receiver's  frequency  measurement  performance. 

The  undeflected  laser  light  is  (the  zero-order  output  from  the  Bragg-cell) 
described  by: 

J =  Lo60  casja'ir£x*b  +■  'Y'toJ  Eq.  11. 

where 

*  undeflected  laser  light  signal. 

L.Cf)  *  amplitude  or  intensity  of  undeflected  laser  light  signal. 

*  laser  light  frequency 

Or‘(fc')  *  phase  variations  on  laser  light  signal. 

The  diffracted  laser  information  signal  is  quantified  mathematically 


Eq.  12. 


Sl,(£)  =. L/t~)  cosj2'ir(?1+¥1')t  +  SltTJ 


where 

jCtl  *  Diffracted  laser  information  signal. 

=  Amplitude  or  intensity  of  diffracted  laser  information  signal. 

4*  =  Frequency  of  information  signal. 

^  (£)  =  Phase  variations  of  information  signal. 

Equation  12  reveals  that  the  diffracted  laser  output  signal  is  frequency 
shifted.  The  frequency  shift  is  produced  by  the  laser  light  signal  and  acoustic 
signal  interaction  in  the  Bragg-celL 

The  Bragg-cell  receiver  is  conceptually  a  parallel,  multi-channel  spectrum 
analyzer,  which  determines  the  Fourier  magnitude  spectrum  of  the  input  IF 
signal.  The  electronic  circuit  model  for  the  Bragg-cell  receiver  is  depicted 
in  Figure  3.  The  IF  signal  input  is  applied  to  a  parallel  narrowband  filter  bank 
with  an  associated  diode  detector  and  Alter  bank.  Each  output  signal  is  a 
power  or  energy  indicator  of  the  IF  signal  spectral  content  in  each  selected 
narrowband  Alter. 

PRACTICAL  BRAGG-CELL  RECEIVERS 

The  ideal  Bragg-cell  receiver  is  characterized  in  the  previous  section, 
with  the  measured  frequency  being  derived  from  Equations  9,  10,  and  12. 
This  section  will  identify  and  briefly  present  some  key  differences  between 
practical  Bragg-cell  receivers  and  the  ideal  bragg-cell  receiver.  The  Bragg-cell 
receiver,  shown  in  Figure  1,  is  a  reference  for  this  discussion. 


BRAGG  CELL 

The  ideal  Bragg-cell  is  assumed  to  have  an  undistorted  acoustic  signal 
replica  (at  acoustic  frequency)  of  the  IF  input  signal  propagating  through  the 
cell.  The  AO  transducer  provides  this  electronic  IF  signal  to  acoustic  signal 
conversion.  The  AO  transducer  has  transfer  function  nonlinearities,  which 
produces  unwanted  and  added  acoustic  signal  components.  Simple  sinusoidal 
electronic  signals  are  subject  to  minor  variations,  while  complex  and  wideband 
signals  may  be  significantly  distorted  by  the  nonlinear  transducer  action. 

The  ideal  Bragg-cell  is  assumed  to  provide  a  perfect  acoustic  signal 
termination,  at  the  far  end  of  the  cell.  In  general,  the  acoustic  signal 
termination  is  not  perfect,  which  cause  low  level  acoustic  signal  reflections 
and  distortion  generated  signal  components  to  propagate  back  through  the 
cell.  These  reflected  acoustic  signals  are  considered  as  additional  distortion 
sources,  hence,  variations  in  the  diffracted  laser  information  signal  can  be 
expected. 

Actual  distortion  effects  are  determined  by  AO  transducer  nonlinearities, 
reflected  acoustic  signal  levels,  total  acoustic  signal  parameter  descriptors 
and  Bragg-cell  properties.  The  complete  acoustic  signal's  instantaneous 
frequency  is  determined  by  amplitudes,  frequencies,  and  phases  of  desired 
acoustic  signal,  undesired  reflected  acoustic  signals  and  nonlinearly  generated 
and  unwanted  acoustic  signals.  Frequency  variations  produced  by  instantaneous 
frequency  changes  will  produce  angle  variations  in  the  diffracted  laser 
information  signal;  thus,  frequency  measurement  errors  in  the  Bragg-cell 
receiver. 
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LASER  SOURCE 

The  laser  light  source  is  ideally  assumed  to  be  constant  amplitude  and 
fixed  frequency  signal.  Actual  laser  sources  can  have  some  minor  frequency 
shifts.  Minor  frequency  shifts  will  results  in  small  angle  variations  in  the 
diffracted  laser  information  signal.  (See  Equation  3.) 

Frequency  variations  will  widen  the  effective  spotsize  and  change  the 
angle  of  the  diffracted  laser  beam,  which  means  the  light  energy  may  be  spread 
over  more  detector  elements  in  the  linear  detector  array. 

The  laser  light  signal  does  not  have  a  constant  intensity  or  amplitude. 
The  amplitude  is  quantified  by  deterministic  and  random  signal  components. 
These  amplitude  variations  will  result  in  diffracted  laser  information  signal 
intensity  level  variations  at  the  detector  plane. 

OPTICS 

Beam  expander  optics  and  diffraction  signal  focusing  optics  are 
fundamental  modules  in  the  Bragg-cell  receiver.  Beam  expander  optics  shape 
the  coherent  laser  beam  to  illuminate  one  entire  side  of  the  Bragy-celL  This 
beam  expansion  provides  a  long  path  interaction  between  the  acoustic  signal 
and  laser  light  within  the  Bragg-celL  Beam  expander  optics  and  laser  source 
are  positioned  so  that  the  laser  light  strikes  the  cell  at  the  Bragg  angle. 
Nonlinearities,  positional  variations,  and  rotational  variations  in  the  beam 
expander  optics  can  produce  amplitude  and  angle  variations  in  the  diffracted 
laser  information  signal.  These  errors  are  usually  very  small  for  well-designed 
optics. 

Focusing  optics  are  used  as  the  so-called  Fourier  transform  lens.  The 
Bragg-cell  is  positioned  in  the  front  focal  plane  of  the  lens.  The  photodetector 


array  is  positioned  at  the  back  focal  plane  of  the  lense.  (Reference  3.)  Again, 
nonlinearities  and  location  variations  in  the  focusing  optics  variations  can 
produce  intensity  and  angle  variations  in  the  diffracted  laser  information  signal. 
These  errors  are  also  very  small  in  well-designed  optics. 

DETECTOR  ARRAY 

The  ideal  detector  array  is  assumed  to  be  a  long  linear  array  with  high 
frequency  resolution  and  sensitivity.  Practical  Bragg-cell  receivers  use  long 
linear  photodiode  arrays,  where  typically  512  or  1024  photodiodes  are  closely 
spaced  along  the  array  length. 

Figure  4  presents  an  illustration  of  two  diffracted  laser  information 
signals  striking  the  photodetector  array  at  two  different  diodes  (two  different 
frequencies.)  In  actual  practice,  the  light  spotsize  for  each  beam  is  larger 
than  the  area  of  one  photodiode;  hence,  the  light  spills  onto  the  area  between 
detectors  and  also  onto  adjacent  detectors. 

Diffracted  laser  information  signal  angle  variations  (frequency  variations) 
produce  energy  spread  over  several  or  many  photodiode  cells.  Each  photodiode 
cell  is  used  to  determine  the  signal  energy  in  resolved  and  calibrated  frequency 
space.  The  photodetector  array  is  electronically  readout  to  determine  frequency 
and  amplitude  information  for  each  diffracted  laser  light  information  signal. 
The  photodetector  array  is  functionally  and  practically  described  as  a  frequency 
sampling  unit,  with  512  to  1024  frequency  bins.  Each  photodiode  determines 
the  elemental  frequency  resolution. 

Bragg-cell  receivers  can  accurately  measure  frequencies  of  time  coincident 
or  time  overlapped  signals,  which  is  a  big  advantage  over  other  receivers, 
such  as  the  IFM  receiver.  These  receivers  are  available  with  bandwidths  to 
approximately  1  GHz,  and  frequency  resolutions  of  100  KHz  to  10  MHz.  The 


8ragg-ceil  receiver's  dynamic  range  is  approximately  30  to  40  db.  -vith  the 
photodetector  array's  dynamic  range  being  the  critical  limiting  factor. 

PROBABILITY  DENSITY  FUNCTION  ANALYSIS 
BRAGG-CELL  RECEIVER 

The  prooaoillty  density  function  (pdf)  of  the  Bragg-cell  receiver's  measured 
frequency  parameter  is  considered  in  this  section.  The  pdf  characteristic 
is  developed  for  a  first  order  Bragg-cell  receiver  model  under  high 
signal-to-noise  ratio  conditions.  The  IF  signal's  frequency  characteristic  is 
assumed  to  have  a  constant  pdf  or  a  Gaussian  pdf. 

Ideal  and  practical  Bragg-cell  receivers  are  discussed  earlier  in  this  report, 
and  characteristics  equations  are  provided  for  the  ideal  receiver.  The  Bragg-cell 
and  AO  transducer  are  key  producers  of  first-order  distortion  effects.  Distortion 
effects  include  AO  transducer  nonlinearities  and  reflected  acoustic  signals 
returning  from  the  far  end  of  the  Bragg-celL  This  analysis  assumed  that  acoustic 
signals  generated  by  AO  transducer  nonlinearities  are  significantly  larger  than 
reflected  acoustic  signals. 

The  ideal  Bragg-ceil  receiver  contains  a  single  acoustic  signal 

Eq.  13. 

which  results  in  the  diffracted  laser  information  signal  given  in  Equation  12. 
fjis  the  frequency  of  the  acoustic  signal,  which  is  directly  related  to  the 
frequency  of  the  input  IF  signaL  (See  Equation  ?.) 

AO  transducer  nonlinearities  are  producing  undesired  acoustic  signals 


in  the  Bragg-cell.  Undesired  or  unwanted  acoustic  signals  are  described  as 
a  cosinusoidal  series  of  N  components  (harmonics  and  intermodulation 
components  at  a  significant  energy  level)  which  are  characterized  as: 


=  U.(t) 


Eq.  14. 


where 

U  (_■£")  =  Total  unwanted  acoustic  signals  in  Bragg-celL  The  AO 
transducer  will  generate  harmonics  of  f^  and  other 
intermodulation  frequency  components. 

=  Amplitude  of  n  th.  unwanted  acoustic  signal  component. 

»  Frequency  of  n  th.  unwanted  acoustic  signal  component. 

Harmonics  of  f1  and  intermodulation  frequency  components 
are  present. 

=  Phase  variations  of  n  th.  unwanted  acoustic  signal  component. 
The  total  acoustic  signal  in  the  Bragg-cell  is  found  by  summing  Equations  13 


and  14: 


^(.t)  +  U-C*) 


Equation  IS  can  be  rewritten,  with  a  lot  of  work,  as: 

U)(f)  a  2(f)  AW  cm  fair  A,*  +  <*,(*)+ /3(f)] 


Eq.  15. 


Eq.  16. 
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The  Instantaneous  frequency  variation,  resulting  from  a  ( t )  phase 
variations  In  Equation  IS,  is 


*  .  '  (»* 
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Eq.  IT. 


The  instantaneous  frequency  variation  is  a  very  complex  function  of 


amplitudes,  frequencies,  and  phases  of  all  unwanted  acoustic  signal  components 


in  the  Bragg-celL  Equation  1?  can  be  rewritten  using  the  Fourier  series 


expansion  as  a  fixed  frequency  plus  a  series  of  harmonically  related  frequency 


components.  A  first-order  approximation  is  the  reduction  of  Equation  17  to 


a  fixed  frequency. 


=  *¥*s 


Eq.  18. 


fg  is  used  to  approximate  the  resultant  frequency  component,  which 
is  produced  by  unwanted  acoustic  signal  components  propagating  in  the 


Bragg-celL 


Resultant  phase  variations  are  approximated  as: 


/SW 


=  2-TT 


Eq.  19. 
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The  total  acoustic  signal  (Eq.  16)  can  be  approximated  as: 


UJ  (■£)  =  -Z(t)  A(t)  cash  ir¥,t  +  +  of ,  (*TJ  m. 


20. 


The  total  instantaneous  frequency  of  w  ( t )  (Equation  20)  is: 


fyt)  =  ¥,  +  % 


Eq.  21. 


Equation  12  can  be  easily  changed  to  approximate  unwanted  acoustic 
signals.  The  modified  Equation  12  results  in  the  following  equation: 

4,Cf)  =  L>)  COSJW&  +  *,  +  $’(*)].  Eq.  22. 

The  actual  displacement  of  the  diffracted  laser  information  signal  (Equation 
22)  relative  to  the  undeflected  laser  light  is 

d  —  - rr - -  •  Eq.  23. 

Equation  23  reveals  that  d  is  a  random  variable  since  it  is  a  summation 
of  two  random  variables  f  1  and  f2>  Conceptually,  the  spot  size  of  the 
diffracted  laser  information  signal  is  increased  by  the  unwanted  acoustic  signal 
components.  The  actual  displacement  d  is  directly  calibrated  to  a  photodiode 


cell  in  the  linear  detector  array,  which  corresponds  to  a  specified  frequency 
bin. 
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1 1  and  f  are  assumed  to  be  incependent  random  variables  for  this 
analysis.  Pdf  characteristics  of  i  1  ard  f  2  are  assumed  to  be  either  constant 
or  Gaussian.  The  pdf  characteristic  of  d  ’  is  the  convolution  of  pdf  (  s  1 ) 
and  pdf  (f^  ) ,  if  the  random  variables  are  independent  (Reference  4.) 

p00 

Pd?(d')  =  pd$(dU,)Pd^(^  df,  E**. 


where 


d'-  cL  7 


Each  frequency  bin  in  the  linear  photodiode  array  exactly  corresponds 
to  Ad  bin  coverage  in  displacement  space. 

The  output  frequency  of  the  detector  array  is 


?D=-  cl 


Substituting  Equation  25  into  Equation  23,  the  result  is: 


Eq.  25. 


Eq.  28. 


The  Bragg-cell  receiver's  measured  output  frequency  can  be  stated  in 
terms  of  the  input  IF  signal's  frequency  by  using  Equations  28  and  a  modified 
version  of  Equation  7.  The  output  frequency  is: 


» 


•T^  =  The  equivalent  IP  frequency  of  unwanted  acoustic  signal 
components  approximated  as  acoustic  frequency  1 2  . 

£  is  a  linear  combination  of  two  random  variables  fi  and  f^2 
,  which  are  assumed  to  be  independent.  The  pdf  (fQ)  is  the  convolution 
of  pdf  (f  )  and  pdf  (  £*  ) ,  which  is  expressed  mathematically  as: 


(?o)  — ^X*)  P^  (£xp)  ^3 


Eq.  28. 


The  pdf  ( f  Q )  characteristic  is  functionally  divided  into  sampled 
bins  by  the  photodiode  cells  in  the  linear  detector  array.  This  frequency  domain 
binning  process  is  equivalent  to  discrete  sampling  of  the  pdf  ( f  ) 

characteristics  in  A£  Q  bin  widths.  A£Q  is  the  basic  frequency 
resolution  of  the  Bragg-cell  receiver  at  the  linear  detector  array.  The  actual 
pdf  characteristic  of  the  Bragg-cell  receiver's  frequency  parameter  is  a  discrete 
characteristic,  which  is  a  function  of  the  pdf  of  the  IF  signal  frequency  and 
pdf  characteristic  of  the  equivalent  unwanted  signal  components. 
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SIGNAL  PROCESSING  ANALYSIS 


INTRODUCTION 

Signal  processing  analysis  results  are  presented  in  this  section.  The 
Bragg-cell  receiver  model  is  implemented  on  a  Hewlett  Packard  Integral 
computer.  Computer  software  program  listings  are  contained  in  Appendix 
A  of  this  report. 

Pdf  analysis  of  the  Bragg-cell  receiver's  frequency  parameter  is 
performed,  assuming  constant  and  Gaussian  pdf  characteristics  for  the  input 
signal's  IF  frequency.  Pdf  analysis  is  performed  for  selected  frequency 
resolutions  of  the  Bragg-cell  receiver.  Also,  statistical  analysis  are  performed 
on  output  pdf  signals  and  summary  results  are  presented. 

BRAGG-CELL  RECEIVER  MODEL  RESULTS 

The  Bragg-cell  receiver  is  math  modeled  in  the  previous  section. 
First-order  model  implementation  is  shown  to  be  a  convolution  of  pdf  (fjp) 
and  pdf  (f *)  for  the  Bragg-cell  receiver's  output  frequency  pdf  characteristic. 
Equation  28  is  the  final  descriptive  equation  for  the  output  frequency  pdf 
characteristic. 

Bragg-cell  receiver  signal  processing  results  are  summarized  in  Table 
1.  Signal  processing  analysis  programs  are  presented  in  Appendix  A.  Figures 
5  through  70  contain  detailed  pdf  plots  for  the  Bragg-cell  receiver's  output 
frequency  parameters.  These  plots  are  made  for  selected  pdf  characteristic 
of  fjp  and  f^  and  selected  Bragg-cell  receiver  frequency  resolutions. 

Table  1  defines  the  input  signal  or  fjp  in  terms  of  DF1  or  DF1  and 
Oj.  The  pdf  (fjy)  is  assumed  to  be  a  constant  if  only  a  value  for  DF1  is  given. 
DF1  (MHz)  indicates  the  frequency  excursion  of  the  input  frequency  (fjp) 

For  a  Gaussian  pdf  (f^.),  DFl(MHz)  and  o^MHz)  is  selected  as 
1.166  MHz  or  4.166  MHz  for  Gaussian  pdfs. 


TABLE  1  CONT'D 
BRAGG-CELL  RECEIVER 
SIGNAL  PROCESSING  RESULTS 
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The  distortion  signal  or  f 2  is  Quantified  in  terms  of  OF2  or  DF2  and 
o 2*  The  pdf(f  j)  is  assumed  to  be  a  constant  if  only  a  value  for  DF2  is  given. 
DF2(MHz)  indicates  the  frequency  excursion  of  the  distortion  signal  f *. 

For  a  Gaussian  pdf(f  2),  DF2(MHz)  and  <?2(MHz)  indicates  six  times 
sigma  and  standard  deviation,  respectively.  DF2  is  selected  as  2,  5,  9,  15, 
and  24  MHz.  a^iMHz)  is  selected  as  0.333,  0.833,  1.5,  2.5  and  4.0  MHz 
for  Gaussian  pdfs. 

Bragg-cell  receiver  frequency  resolution  (MHz)  is  chosen  to  be  0.5, 
1.0,  2.0  or  5.0  MHz.  These  frequency  resolution  selections  are  readily 
expected  with  current  Bragg-cell  receiver  technology. 

The  last  three  columns  of  Table  1  contain  statistical  results  of  the 
Bragg-cell  receiver's  output  frequency.  Sigma  (MHz),  skew  and  kurtosis 
are  computed  for  each  output  frequency  pdf  characteristic  shown  in  Figure 
5  through  70. 

The  pdfif J  characteristic  is  a  staircased  trapezoidal  function  for  pdfffjp) 
=  constant  and  pdfif^)  =  constant.  Figures  5  through  22  show  numerous 
examples  of  pdfif^)  characteristics  for  selected  input  signal  and  distortion 
signal  constant  pdf  characteristics.  Pdfif  ^  is  observed  to  significantly  depart 
from  a  constant  characteristic  and  approach  a  staircased  triangular 
characteristic  as  OF2  approaches  DF1. 

Figures  23  through  34  reveal  output  frequency  pdfif^)  characteristics 
with  input  frequency  characteristics  assumed  to  be  Gaussian  pdf  and  distortion 
signal  frequency  characteristics  assumed  to  be  constant  pdf.  These  pdf 
characteristics  are  staircased  Gaussian  characteristics.  The  pdf  characteristic 
departs  from  the  Gaussian  characteristics  as  the  Bragg-cell  receiver's 
frequency  resolution  is  decreased. 


Figures  23  through  34  reveal  output  frequency  pdf(f^)  characteristics 
with  input  frequency  characteristics  assumed  to  be  Gaussian  pdf  and  distortion 
signal  frequency  characteristics  assumed  to  be  constant  pdf.  These  pdf 
characteristics  are  staircased  Gaussian  characteristics.  The  pdf  characteristic 
departs  from  the  Gaussian  characteristics  as  the  Bragg-cell  receiver's 
frequency  resolution  is  decreased. 

Figures  35  through  52  reveal  output  frequency  pdf(f^)  characteristics 
with  input  frequency  characteristics  assumed  to  be  constant  pdf  and  distortion 
signal  frequency  characteristics  assumed  to  be  Gaussian  pdf.  For  most 
example  plots  the  long  constant  pdf  characteristic  of  the  input  frequency 
tends  to  dominate  the  overall  pdf  characteristic.  The  pdf  characteristics 
are  approximately  a  staircased  trapezoidal  characteristic. 

The  pdf(f^)  characteristic  is  a  staircased  Gaussian  function  for  pdf(fjp) 
=  Gaussian  and  pdf(f  *)  a  Gaussian.  Figures  53  through  70  show  numerous 
examples  of  pdf(f^)  for  selected  input  signal  and  distortion  signal  frequency 
characteristics.  Pdfif^)  is  observed  to  be  close  to  a  Gaussian  characteristic 
for  all  DF2  and  DF1  selections.  For  low  frequency  resolution  in  the 
Bragg-cell  receiver,  the  pdf  characteristic  departs  from  a  Gaussian 


characteristic. 
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CONCLUSIONS 


Bragg-cell  receivers  can  accurately  measure  frequencies  of  time 
coincident  and  time  overlapped  RF  signals,  which  is  a  big  advantage  over 
other  receiver  types.  Bragg-cell  receivers  are  available  with  bandwidths 
to  approximately  1  GHz,  and  frequency  resolutions  of  100  KHz  to  10  MHz. 

The  Bragg-cell  receiver's  frequency  parameter  is  analyzed  and  results 
are  presented  in  this  report.  Probability  density  function  analysis  and 
statistical  analysis  results  are  presented  for  this  receiver  for  selected 
frequency  resolution  capabilities. 

The  first  order  math  model  is  developed  for  the  Bragg-cell  receiver. 

The  input  signal's  frequency  is  assumed  to  be  a  constant  pdf  or  Gaussian 
pdf  characteristic.  The  primary  distortion  signal  is  generated  by  nonlinear 
acoustic  transducer  characteristics.  The  distortion  signal  is  also  assumed 
to  have  a  constant  pdf  or  Gaussian  pdf  characteristic. 

Staircased  pdf  characteristics  are  observed  for  all  pdf  frequency 
characteristics.  The  staircased  function  is  introduced  by  the  frequency 
sampling  of  the  photodiode  detector  array  used  in  the  Bragg-cell  receiver. 

The  output  frequency  pdf  characteristic  is  observed  to  depart  from  the  input 
signal's  frequency  pdf  characteristic  for  many  of  the  wideband  RF  signal 
cases  and  wideband  distortion  signals  produced  by  the  AO  transducer. 

Distortion  signals  with  large  DF2  values  can  significantly  change  the 

output  frequency's  pdf  characteristic.  Also,  the  frequency  sampling 

(or  detector  frequency  resolution)  at  the  detector  can  further  distort 

the  pdf  characteristic.  Low  frequency  resolution  at  the  detector  produces 

the  most  significant  distortions.  EW  emitter  classification  and  EW 

direction  finding  systems  can  be  affected  by  Bragg-cell  receiver  distortions. 

These  distortions  are  especially  important  for  wideband  radar  emitters. 


ViVi'jr  Ji 


u*  »-•  >U’4. 


Math  model  validation  is  recommended  for  the  Bragg-cell  receiver. 
Model  validations  are  readily  performed  using  an  actual  Bragg-cell  receiver 
and  associated  test  equipment. 
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APPENDIX  A 
BRAGG-CELL  RECEIVER 
CONVOLUTION  PROGRAM 


1 0  !  “BRAGG'*  PROGRAM 

23  '  BRAGS- CELL  RECEIVER  ANALYSIS 

33  !  CONVOLUTION  OF  F  5  ARRAY  ANO  F2  ARRAY 

43  1  THIS  VERSION  PERMITS  STORAGE  OF  CONVOLVED  ARRAY  UITH  ITS  FREQUENCY  AXIS 

S3  !  STORAGE  TAKES  PLACE  PRIOR  TO  PLOTTING— A  STRING  NAMED  INFOS  IS  ALSO 

63  !  STORED— SHOULD  CONTAIN  MINIMUM  INFORMATION  ON  STORED  DATA 

73  I  PLOT  BOTH  INPUT  ARRAYS  ON  K  AXIS 

30  !  PLOT  CONVOLUTION  ON  FREQUENCY  AXIS 

33  !  LABEL  TOTAL  VALUES  FOUNO  IN  CONVOLUTION  ARRAY 

133  !  ARRAYS  MAY  BE  ANY  COMBINATION  OF  FLAT/GAUSSIAN 

113  :  CECSMSER  1388 

123  'ISP  “CO  YOU  WANT  TO  SAVE  CONVOLVED  ARRAYS?**  3  LINPUT  “ENTER  Y/N",XPS 
133  ! IF  KPSC 1 . 1 3*“Y“  THEN  30SUB  KEEP 

143  CIS?  “ENTER  INTEGER  LENGTH  OF  LONG  ARRAY"  8  INPUT  DF1 
1S3  IDISP  "ENTER  LENGTH  OF  SHORT  ARRAY"  3  INPUT  0F2 
183  DISP  “MAKE  THESE  ENTRIES  IN  MHZ" 

173  01 3P  “ENTER  DELTA  FREQUENCY  FOR  THE  LONG  ARRAY"  8  INPUT  DFFl 
183  0I3P  “ENTER  DELTA  FREQUENCY  FOR  THE  SHORT  ARRAY"  8  INPUT  OFFZ 
133  0F2»CF1/<0Fr 1/0FF2)  8  N-CF1*Z  8  3GF1-0FF 1*1333330/8  8  SGFLO-QFFZ* 1303033/8 
233  OISP  “OF* »"  i  DF1 1 “  DF2-  *  lDF2  i '*  N»"jN  8  NMR1-I /<  SGF1  *SQR<  2*PI )  > 

213  OPTION  BASE  1  8  RAO  8  NMRLG-  1  /<  SGFLQ*SCR(  2*PI  > )  8  MU-533033333 

223  DIM  PX<  SCO ) . PH<  533 > . P Y( 499 ) . TSt  32 1  ,  T 1  St  52 3  . FX( 533 > . Xt  533 > , Y<  533 > . INFOSt  831 

233  LINPUT  “ENTER  -<83  CHARACTERS  DESCRIBING  DATA". INFOS 

243  REDIM  PYt N- 1 ) ,FX< N-1 ) ,PX( N) ,PH< N) 

253  FOR  K-l  TO  N  !  FILL  INPUT  ARRAYS  UITH  ZEROES 
253  PX( X  )*3  8  PH<  X ) -3 
273  NEXT  K 

233  OISP  “  *****PICX  COMBINATION*-***" 

293  OISP  "  J-BOTH  ARRAYS  FLAT" 

303  OISP  "  2-SHORT  FLAT,  LONG  GAUSSIAN" 

313  OISP  "  3-LONG  FLAT.  SHORT  GAUSSIAN" 

323  OISP  "  4- BOTH  GAUSSIAN" 

333  OISP  “ENTER  YOUR  CHOICE"  8  INPUT  PIC 
343  IF  PIC-2  THEN  SHRT 
353  IF  PIC-4  THEN  60SU8  AR2 
363  LSTEP-1 

373  FOR  K«1  TO  M\2  I  MAKES  FLAT  LONG  ARRAY 
383  PH(K)«1/(0F1-1 ) 

393  NEXT  X 

433  SHRT:  I  MAKES  FLAT  SHORT  ARRAY 

413  IF  PIC-3  THEN  30SU8  AR3 
423  SSTEP-1 

433  FOR  K-N/4-CF2/Z  TO  N/4+CF2/2 
443  PX(K)»1/(DF2-1 ) 

4S3  NEXT  K 

463  IF  PIC-2  THEN  30SUB  AR2 
473  CNVLV: 

483  YTOT-3 

493  FOR  K-l  TO  N-1 
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so®  vsum«o 

S10  STL-SO0-CFF1  a  STPL-S0O+OFF1  a  STZL*< STPL-STL)/(  N- 1 )  3  IVAL-STL+ST2L/2 

52®  FOR  J-t  TO  K 

S3®  YSUM-YSUM+PXC  J)*PH(X-J+1 ) 

540  PY<  X  >«YSUM/<  STZL*  I0OO00®) 

SS®  NEXT  J 

55®  IF  X  -TOO  SO-0  THEN  OISP  "X-"^ 

5*7©  M£aT  X 

S3®  FOR  X«1  TO  M-!  !  TAKE  FREQ  ARRAY  FGR  XAXIS 

S3®  FX1X1-IVAL 

50®  IVAL-IVAL+STZL. 

SI®  NEXT  X 
S20  SUMY*® 

S3®  FOR  X*1  TO  N-2  'INTEGRATE 
64®  XS*<FX(X+1 )-FX<X 5 >*100000® 

SS®  YS-<  PY{ X  >+PY<  X+1 > )/2 
SS®  3UMY*XS*YS+SUMY 
S7®  NEXT  X 

630  01 SP  "SUM--; SUMY 

S3®  ! IF  XP*t 1,11 -"Y"  THEN  30SU8  XP1 


70®  CHOICE:  !  PICK  PLOT 

71®  OISP  *  1-PLOT  BOTH  ORIGINAL  ARRAYS  SAME  AXIS- 

72®  01 SP  -  2-PLOT  CONVOLUTION  ON  FREQUENCY  AXIS- 

73®  OISP  “  3-EXIT  PROGRAM- 

740  OISP  “  4- STORE  OATA“ 

750  OISP  “  ENTER  YOUR  CHOICE"  Q  INPUT  CH 

7S3  IF  CH-4  THEN  GOSUB  XEEP 

770  IF  CH«1  THEN  73® 

73®  LINPUT  -EXTERNAL  PLOTTER?  Y/N“,HP* 

73®  IF  CH-3  THEN  BEEP  «  OISP  "OONE"  0  STOP 
30®  IF  CH-1  THEN  GOSUB  PLT1 
81®  IF  CH-2  THEN  GOSUB  PLT3 

82®  PLOTTER  IS  1  0  PEN  -1  B  GCLEAR  B  LOCATE  30,170.25,3® 

33®  Y1-®  B  Y2-AMAX1 Y)*  1 . 1  B  Xl-STL  0  X2-STPL 

84®  START: 

35®  SCALE  X1.X2.Y1.Y2 

86®  FXO  2  B  LAXES  ( X2-X1 )/1®,< Y2-Y1 )/10,X!  ,Y1 
37®  XAXIS  Y2 ,<  X2-X1 )/10,X1 ,X2 
88®  YAXIS  X2,<Y2-Y1 J/10.Y1 ,Y2 

83®  MORCAT:  I  ADO  MORE  OATA  SAME  SCALE 

30®  IF  CH-1  THEN  35® 

31®  IF  CH-4  THEN  35® 

32®  OISP  "ENTER  8RAGG- RECEIVER  FREQUENCY  RESOLUTION- 
33®  INPUT  RES 

34®  KR£S-INT(  RES/1 Z»QFF1/N) ) 

350  FOR  X-1  TC  TOP 

36®  PLOT  X( K ) ,Y<  X )  B  IF  CH-1  THEN  104® 

370  IF  CH-4  THEN  1040 
38®  FOR  L-1  TO  XRES 
33®  IF  K+L>«TCP  THEN  1050 
1000  PLOT  X(X+L),Y(K) 

101®  NEXT  L 
102®  K-K+L-2 
103®  L-1 
1040  NEXT  K 

105®  IF  HP*-"Y-  THEN  EPLOT 

1060  IF  CH-1  THEN  S0SU8  PLT2 

1070  IF  CH-4  THEN  PLOUT  103 
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5080  IMAGE  “CONSTANT  PDF  FOR  FT.  OF t .00. 3. -  CONSTANT  P OF  FOR  FZ  .  CFFZ»“.00.0 

1090  IMA6E  “CONSTANT  POF  FOR  FT.  OFT »“ ,00.0.“  GAUSSIAN  PDF  FOR  FZ  .  CFFZ»“.DD.D 
*• 

5  100  IMAGE  “GAUSSIAN  POF  FOR  FI.  0F1-“  .00.0.“  CONSTANT  POF  FOR  FZ  .  GFF2-“,D0. 

0.“ 

15  50  IMAGE  “GAUSSIAN  PDF  FOR  FI.  CF5-“  ,00.0,"  GAUSSIAN  POF  FOR  FZ  .  OFFZ-". 00. 
0." 

1TZ0  7*«“3RAGG-CSLL  RECEIVER  POF  ANALYSIS"  8  CSIZE  3.3, T 

IT  30  LORG  5  8  MOVE  XZ-<  X2-X1  )/Z  ,YZ+<  Y2-Y1 >/3  8  LABEL  USIN6  “X"  |  T* 

1 540  CSIZE  3 

11S0  MOVE  X2-<  X2-X1 >/2 , Y2+<  Y2-Y1 )  / 1 3  8  IF  PIC-1  THEN  LABEL  USIN6  5080  «  OFF 5 .OFF 
Z  8  CSIZE  3 

5  ISO  IF  PIC-3  THEN  LABEL  USING  5090  t  OFF? ,DFF2  8  CSIZE  4 

5570  IF  PIC-Z  THEN  LABEL  USING  5100  j  0FF1 .OFFZ  8  CSIZE  4 

1180  IF  PIC-4  THEN  LABEL  USING  1 1 50  »  OFF!  , OFFZ  8  CSIZE  4 

1 1 90  MOVE  X2-<  X2-X5 )/Z.Y1-< Y2-Y1 >/S  8  LABEL  USING  “X"  »  “FO  IN  MHZ" 

1200  MOVE  XI -< X2-X1 )/7, Y2-< Y2-Y1 )/2  3  CES  8  LOIR  30  8  LABEL  USING  “X"  j  "Pdf(FO) 
X  IEB" 

1250  IMAGE  “TOTAL-". 00. CCC 

1220  MOVE  XI ,Y1-<YI-Y1)/S  8  LOIR  0  8  LORG  2  8  LABEL  USING  52tO  j  SUMY 
1230  PLCUT: 

5240  L INPUT  “GUMP  GRAPHICS?  Y/N“,DG* 

12S0  IF  CG*-"Y"  THEN  DUMP  GRAPHICS 
12S0  EPLOT: 

1270  HP*»“N“ 

1230  LINPUT  “PLOT  ON  EXTERNAL  DEVICE?  Y/N" ,PL£» 

1230  IF  PLE*-“Y"  THEN  PLOTTER  IS  70S  8  PEN  1  8  LOCATE  30.110.29.89  8  SOTO  START 
1300  IF  CH-4  OR  CH-Z  THEN  CHOICE 
1370  BEE?  150.300  8  OISP  “DONE" 

1 320  END 

1330  PLT1 :  1  PLOT  SHORT  ARRAY 

1340  REDIH  X(N).YtN) 

1350  IF  CH-1  THEN  NULT-t  ELSE  MULT- 1 00000000 
1360  FOR  K«1  TO  M  I  FILL  X  AND  Y  ARRAYS 
1370  XIX )«K  8  Y(K)-PX(K)-MULT 
1380  NEXT  X 

1390  X1-1  8  X2-N  8  Y1-0  8  Y2-AMAX( Y1-1 . 1  8  TOP-N  8  LINE  TYPE  1 
1400  PLOTTER  IS  1  8  PEN  -1  8  SCLEAR  8  LOCATE  30.170,25,80 
1410  GOTO  START 
1420  RETURN 

1430  PLTZ:  1  PLOT  L0N6  ARRAY  SAME  SCALE 

1440  T*-“ INPUT  ARRAYS" 

1450  IMAGE  "CONSTANT  ARRAY  FOR  FI,  0F1-" ,00.0,"  CONSTANT  ARRAY  FOR  F2  .  0FF2-". 
00.00 

1480  IMAGE  "CONSTANT  ARRAY  FOR  FI.  OFl-“,DO.D,"  GAUSSIAN  ARRAY  FOR  FZ  .  OFFZ-", 
00.00 

1470  IMAGE  "GAUSSIAN  ARRAY  FOR  FI,  DFt-“.DC.O,“  CONSTANT  ARRAY  FOR  F2  ,  0FF2-". 
00.00 

1480  IMAGE  "GAUSSIAN  ARRAY  FOR  FI.  DF!«“.DO.D,"  GAUSSIAN  ARRAY  FOR  FZ  .  OFFZ-", 
00.00 

1490  LORG  5 

1500  MOVE  Xr-<X2-Xn/2.Y2*<  Y2-YU/20  8  IF  PIC-2  THEN  LABEL  USING  1470  |  0FF1.0FF 
2  8  CSI2E  4 

1510  MOVE  X2-( X2-X1 )/Z . YZ+< Y2-Y1 )/Z0  8  IF  PIC-1  THEN  LABEL  USING  1450  I  OFFI.OFF 
2  8  CSI2E  4 

1520  MOVE  X2-( X2-X1 1/2 ,Y2+< Y2-Y1  )/20  8  IF  PIC-3  THEN  LABEL  USING  1480  »  OFFI.OFF 
2  8  CSIZE  4 
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1530  MOVE  XZ-<X2-X1  >/2 .  Y2+<  YZ-YD/ZO  3  IF  PIC-4  THEN  LABEL  USING  1480  i  OFF  1  .OFF 
2  3  CSI2S  4 

1540  LORG  5  3  HOVE  X2-<X2-X1)/2.Y2+<Y2-Y1)/9  3  LABEL  USING  “X“  i  T# 

1550  IF  CH-I  THEN  MULT- I  ELSE  J1ULT- 1 000000000 

1SS0  .10VE  XI -<  X2-X1  1/7, Y2-<  Y2-Y1  )/Z  3  DEG  3  LOIR  30  3  LABEL  USIN6  “K“  i  “VALUE  X 
“;J1ULT 

1570  FOR  X-1  TO  N 
1580  Y<  X )  «PH<  X )  *MULT 
1530  NEXT  X 

1600  CH-4  3  LINE  TYPE  4  3  TOP-N 
1S10  SOTO  HOROAT 
1620  RETURN 

1630  PLT3:  !  FILL  CONVOLUTION  ARRAY 

1640  REOin  X(  N- 1 ) ,  Y<  N- 1 ) 

1650  NULT* 1000000000 

1660  FOR  X-1  TO  N-1 

1670  X(X)-FX<X 5  3  Yt K  >-PY( X )*(1ULT 

1630  NEXT  X 

1630  TOP-N- 1 

1700  RETURN 

1710  AR2:  ! SU8-SENERATE  GAUSSIAN  FOR  THE  LONG  ARRAY 

1720  Stt«flU-3*SGFl  3  3tp«f1U+3*S6Ft  Q  ST2-< Stp-Stt )/< N\2 )  8  IVL-Stt 

1730  FOR  X-1  TO  N\2 

1740  PH<  X )«NMR1 »£XP< ( -< IVL-HU)“2 )/( 2-S6F1 *Z> >  8  PH< K )-PH( X )*STZ 
1750  IVL«IVL>ST2 
1760  NEXT  X 

1770  IF  PIC-4  THEN  AR3 
1780  IF  PIC-2  THEN  CNVLV 
1730  RETURN 

1800  AR3:  ! SU8-GENERATE  GAUSSIAN  FOR  THE  SHORT  ARRAY 

1810  Stt-MU- 3-SGFLO  8  Stp-MU+3*SGFL0  •  STZ»< Stp-Stt >/CN/4+0FZ/Z-< N/4-0F2/Z) >  8  I 
VL-Stt 

1820  FOR  X-N/4-QF2/2  TO  N/4+CFZ/Z 

1830  PX<  X )*NMRLQ*EXPt (t< IVL-MU)“Z)/< Z*SGFLQ“2> )  8  PX(K)-PX(X)-ST2 
1840  IVL-IVL+STZ 
1850  NEXT  X 

>860  IF  PIC-3  OR  PIC-4  THEN  CNVLV 
1870  RETURN 

1880  KEEP:  I  FILENAME  ENTRY 

1830  OISP  “FILENAME  MUST  BE  ENTERED  IN  UNIX  PATH  FORM  /DXSCNAME/FILENAME“ 

1300  L INPUT  “ENTER  DESCRIPTIVE  FILENAME" ,F* 

1310  CREATE  F8. 1 .299*2+33 

1920  ASSI6N#  I  TO  Ft 

1930  PRINT#  1  i  N,INFQt,FX( > ,PY<  ) 

1940  ASSIGN#  1  TO  • 

1950  GOTO  CHOICE 
1960  RETURN 


APPENDIX  A 
8RACC-CELL  RECEIVER 

PDF  GENERATION  &  STATISTICAL  ANALYSIS  PROGRAM 


S  ! 8RAGGP0F  THIS  PROGRAM  COMPOTES  3RAG6-CELL  RCVR.  POF  A NO  STATISTICS 
;3  i  'HIS  PROGRAM  SETS  DATA  FROM  A  CONVOLUOLUTICN  FILS-00E5  CEN  MOMS 

03  'CATA  IS  THEN  PLOTTED  WITH  LABELLED  CENTRAL  MOMENTS.  INT .  ETC - 

33  OPTION  3ASE  : 

43  DIM  FX<  583  * ,PH<  S38 ) , INFO*t  803 ,T*£75] 

53  3 ISP  'PLEASE  USE  FILENAMES  WITH  0  DIGITS  FOR  CF1* 

53  L INPUT  -FILENAME?* ,F* 

73  ASSIGNS  ;  TO  =* 

33  READS  1  i  N  3  N— N— * 

33  SEDIM  FX(  N  !  ,?H<  N ) 

133  READS  !  i  INFO*,FX<  >  ,PH<  ) 

M3  "OR  :<■'  TO  N 
t23  PX<K  >»*XCX  :+>33 
i 33  NEXT  K 
I  A3  ASSIGN*  5  TO  * 

;S3  OISP  3  3ISP  INFO*  3  0I3P 

t so  =!3«RE';s<Fa;  a  p*«ci*u  ,n 

173  IF  PS«-A“  THEN  PIC- I 

133  IF  P*«-3*  THEN  PIC-0 

138  IF  P9-*C*  THEN  PIC-3 

033  IF  P*«*3*  THEN  PIC-+ 

013  0F1»'/AL(F3I0.4I  >  9  3FLC»UAL< F*CS  ,?I > 

023  OISP  “ENTER  RESOLUTION  IN  MHZ*  9  INPUT  RES 
233  XST»INTv'RE3/«'2*0Fl/N) ) 

0*3  CFM-AMAX<FX>-AMIN<FX>  9  FSTPM-QFM/1 N-U-XST  9  OF-FSTPM* 1333388  - 9  SUMZ-3 
ZS3  SUMY-3  9  3UMS-3  9  SUM3-3  9  SUM4-0 
2S3  FOR  K— I  TO  N-1 

273  XS-(FX<K+1  )-PX(K  ))• 1380888  9  YS-< PH<K >+PH<K+1  > )/2 
033  3UMZ-SUMZ+XS-YS 
033  NEXT  K 

333  FOR  K-1  TO  N  STEP  XST 

313  FX  < K )»FX  <  K  > • ! 383333 

303  PH<  X  >-PH<  K //SUMZ 

333  SUMY-SUMY+PH<  K  )«FX.(  K  )  *0F 

3*0  NEXT  K 

353  XMEAN-SUMY 

363  FOR  K-1  TO  M  STEP  XST 

373  SUMS-SUMS+<FX(K  J-XMEAN  )“2*PH< K  1-OF 

380  NEXT  K 

393  UNO -SUMS  9  3IGMA-SCR< VNC > 

433  FOR  K-1  TO  N  STEP  XST 
418  SUM3»SUM3+<  FX<  K  1-XMEAN  >“3*PH<  K  >«0F 
423  SUM4-SUM4+(FX(K l-XMEAN  )*4*PH<K  )*0F 
433  NEXT  K 

448  CM3-SUM3  9  CM4-SUM4 

4S8  SKEU-CM3/UNC* 1 .5  9  KURT-CM4/VNC*2 

463  PLOTTER  IS  1  9  PEN  -1  9  GCLEAR  9  LOCATE  38.178,25.33 

473  Y1-3  9  Y0-AMAX1 PH  5*1 38888888 • I . 1  9  Xt-FX< 1)/ 1930308  9  XO-XS-N/ 1000000/2+530 
480  START: 


r,i  ■■«3  "ll  . 


gL«.k  >,t 


433  SCALE  Xt  .:<:  ,Y1  ,Y2 

533  FXO  3.2  3  LAXE5  tX2-<1  >/10((Y2-Y1  ;/l3,Xt  ,Y1 
5 i 3  XAX is  yz  ,tX2-<1  t/13  ,XJ  ,X2 
523  -'AXIS  X2  YZ-/1  >/10.Y1  ,Y2 

523  MCROAT:  '  AOO  MORE  CATA  SAME  SCALE 

543  FOR  K*l  TC  N  STEP  X3T 
553  IF  K+X3T>N  THEN  533 

SS3  PLOT  FX<K  333333  ,PH<K  )♦  1 33303333  9  PLOT  FXtK+XST  >/ 1333830  ,?H(K  >•  183303880 
573  NEXT  K 

5=3  IMAGE  “CONSTANT  POF  FOR  Ft,  0F1 ■*  .000.0 .*  CONSTANT  PCF  FOR  F2 .  0FF2-* ,000.0 
3 

530  IMAGE  “CONSTANT  PDF  FOR  Ft.  3F1 ■“  .300.3 ,“  3AUSSIAN  POF  FOR  F2  .  CFF2-*  ,000.3 
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533  IMAGE  “3AUS3IAN  POF  FOR  Ft  .  OFt «*  ,300.3 CONSTANT  POF  FOR  F2  .  0FF2-*  ODD. 
30 

St  3  IMAGE  “3AU33IAN  POF  FOR  Ft,  OF! -* .300.3 .*  3AUS3IAN  POF  FOR  F2 ,  DFF2-*  300. 
30  ‘ 

528  Ti-'3RAGG-CSLL  RECEIVER  POF  ANALYSIS*  9  CSI2E  3.3.1 

523  LORG  5  3  MOVE  X2-<  X2-X t  >/2  ,Y2+<  Y2-Y1  )/9  9  LABEL  USIN6  “K*  »  T* 

•O—v  wdii  C 

SS3  MOVE  X2-<X2-Xt  )/2 ,Y2+( Y2-Y1 >/t8  3  IF  PIC-t  THEN  LABEL  USING  550  i  OFt  3FL0  3 
-  C3IZE  4 

553  IF  PIC-3  THEN  LABEL  USING  530  i  OFt  ,3FL0  3  CSIZE  4 

573  IF  PIC-2  THEN  LABEL  USING  530  »  OFt  ,QFLO  3  CSIZE  4 

530  IF  PIC-4  THEN  LABEL  USING  513  «  0F1  .0FL3  9  CSIZE  4 

530  MOVE  X2-<X2-Xt  )/2.Y1-(Y2-Y1  )/S  3  LORG  5  3  LABEL  USING  *K*  i  *Fd  IN  MHZ- 
730  L0R6  5  9  CSIZE  4  9  LORG  S 

713  MOVE  X1-<X2-X1 )/7 ,Y2-< Y2-Y1  )/2  3  0E5  3  LOIR  30  3  LABEL  USING  *K*  i  *Pdf(Fd> 

X  1E3* 

723  LOIR  3  3  RAC 
733  IMAGE  "INTEGRAL-* ,0.000 
743  IMAGE  -XMEAN-* .00D00.000 
753  IMA6E  ‘SIGMA-* ,000.000 
753  IMAGE  ‘SKEW-*  .SCO. 000 
770  IMAGE  ‘KURTOS IS- *,00.000 
780  CSIZE  Z.S  9  LORG  2 

790  MOVE  X1+<X2-X1 )/2S ,Y2-( Y2-Y1  >/ 1 0  9  LABEL  USING  730  «  SUMZ 

800  LABEL  USIN6  740  t  XMEAN/ 1 000000 

810  LABEL  USING  753  I  SIGMA/1 000000 

823  LABEL  USING  763  i  SKEW 

833  LABEL  USIN6  770  t  KURT 

843  LABEL  USING  *K*  i  “RESOLUTION-* |RE3 1  *  MHZ* 

853  LINPUT  “OUMP  8RAPHICS?  Y/N*  ,80M8 
863  IF  80M8I1 ,1 3-*Y*  THEN  OUMP  GRAPHICS 
873  IF  EPL*-*Y*  THEN  900 

880  LINPUT  'PLOT  ON  EXTERNAL  0EVICE7  Y/N* ,EPL* 

890  IF  EPL*«*Y*  THEN  PLOTTER  IS  70S  •  PEN  1  9  LOCATE  30,110,29,39  0  GOTO  START 
903  BEEP  303,433  9  OISP  *CCNE* 
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